Magnetization plateau in S = 3/2 antiferromagnetic Heisenberg chain 

with anisotropy 



On 
On 



O 
o 



Toru Sakai 1 and Minora Takahashi 2 
1 Faculty of Science, Himeji Institute of Technology, Kamigori, Ako-gun, Hyogo 678-12, Japan 
2 Institute for Solid State Physics, University of Tokyo, Roppongi, Minato-ku, Tokyo 106, Japan 

(October 97) 

The magnetization process of the 5=3/2 antiferromagnetic Heisenberg chain with the single-ion 
anisotropy D at T — is investigated by the exact diagonalization of finite clusters and finite-size 
scaling analyses. It is found that a magnetization plateau appears at m — 1/2 for D > D c = 
0.93 ±0.01. The phase transition with respect to D at D c is revealed to be the Kosterlitz-Thouless- 
type. The magnetization curve of the infinite system is also presented for some values of D. 
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One-dimensional antiferromagnets have various quan- 
tum effects observed even in macroscopic measurements. 
The Haldane gap JlJ , which is the lowest excitation gap of 
the ID Heisenberg antiferromagnets with integer 5, was 
also detected as a transition from a non-magnetic state to 
magnetic one in high-field magnetization measurements 
of Ni(C 2 H 8 N 2 )2N02(C104), abbreviated NENP, which 
is an 5 = 1 quasi- ID antiferromagnet. m Recently Os- 
hikawa,Yamanaka and Affleck suggested that even for 
the ID 5 = 3/2 (half-odd integer) antiferromagnet an 
energy gap is possibly induced by a magnetic field and a 
magnetization plateau appears at m = 1/2, which corre- 
sponds to 1/3 of the saturation moment. Their argument 
is based on the analogy to the quantum Hall effect and 
the valence bond solid picture for 5 — 1. || The magneti- 
zation plateau is also predicted in some alternating spin 
chains J|||, but the mechanism depends on the struc- 
ture of the unit cell and the argument for them is not 
necessarily valid for uniform chains. 

For the anisotropic 5 = 3/2 antiferromagnetic chain, 
a variational approach j^] gave the phase diagram of the 
nonmagnetic ground state, while few works were done 
on the magnetic state. However, it is easy to under- 
stand that it should have a magnetization plateau at 
least when the system has the positive and infinitely large 
single-ion anisotropy DJ2j(Sj) 2 - Because in the limit 
(D — > oo ) every site has S? = 1/2 for the ground state 
at m = 1/2 and any magnetic excitations changing it 
into Sj = 3/2 at a site have a gap proportional to D. 
For finite D, however, there is no rigorous proof on the 
existence of the gap at m — 1/2, in contrast to the case 
of m 7^ 1/2 in which the system is proved to be gapless 
by the Lieb-Schultz-Mattis theorem, [^[jj Thus some nu- 
merical tests are important to check the existence of the 
gap and magnetization plateau at m — 1/2. The density 
matrix renormalization group approach ||] revealed that 
the isotropic 5 = 3/2 antiferromagnetic chain is gapless 
even at m = 1/2 and a critical value D c should exist as 
a boundary between the gapless and massive phases. 

In this paper, using the exact diagonalization of finite 



clusters up to the system size L = 14 and finite-size scal- 
ing analyses, we investigate the 5 = 3/2 antiferromag- 
netic Heisenberg chain with the single-ion anisotropy and 
estimate the critical value D c at m = 1/2 and determine 
the universality class of the phase transition with respect 
to D. In addition we present the ground-state magne- 
tization curve extrapolated to the thermodynamic limit 
for some typical values of D. 

Consider the ID 5 = 3/2 antiferromagnetic Heisenberg 
Hamiltonian with the single-ion anisotropy in a magnetic 
field 



H = Hq + Hz, 



(1) 



n 2 



under the periodic boundary condition. For L-site sys- 
tems, the lowest energy of Ho in the subspace where 
y\ = M (the macroscopic magnetization is m — 
M/L) is denoted as E(L, M). Using Lanczos' algorithm, 
we calculated E(L, M) (M = 0, 1, 2, • • • , 3X/2) for even- 
site systems up to L = 14. For finite systems described 
by the total Hamiltonian H, the energy gap of the mag- 
netic excitation changing the value of M by ±1 is given 
by 



A± = E(L, M ± 1) - E(L, M) =F H. 



(2) 



If the system is gapless in the thermodynamic limit, the 
conformal field theory(CFT) gives the asymptotic form 
of the size dependence of the gap as A± ~ 0(1/ L) with 
fixed m = M / L. If we define H + and i/_ as 



E(L, M + 1) - E(L, M) -> H+ (L 
E(L, M) - E(L, M — 1) — » H- (L 



(3) 



H + and H— has the same value and it gives the magnetic 
field H for the magnetization m in the thermodynamic 
limit. On the other hand, if the system has a finite gap 



1 



even in the limit, neither A + nor A_ vanishes for L — > oo. 
It implies that H+ and H- are different. As a result, a 
plateau appears for H- < H < H + at m — M/L in the 
ground-state magnetization curve. 

Since A± includes an undecided parameter H in the 
form (|J), we take the sum A = A + + A_ for the order 
parameter of the finite-size scaling, to test the existence 
of the plateau at m = 1/2. (In the massive case, the gap 
A leads to the length of the plateau in the magnetization 
curve in the thermodynamic limit.) The scaled gap LA of 
finite systems (L = 6 ~ 14) at m — 1/2 is plotted versus 
D in Fig. |l|. For D > 2 the scaled gap obviously increases 
with increasing L, which means that a finite gap exists in 
the thermodynamic limit. For small D around the region 
< D < 1, the scaled gap looks almost independent of 
L. It implies that the system is gapless at a finite region. 
At least the form A ~ 1/L is valid for < D < 0.8 with 
the relative error less than 0.3% for each point. Our 
precise analysis, however, indicates that the LA curves 
for L, and L + 2 have only one intersection in the region 
< D < 2 for each L. Thus the critical point D c can 
be estimated by the phenomenological renormalization 
group equation p2[ 



(L + 2)A L+2 {D')=LA L {D). 



(4) 



We define D c l,l+2 as the L-dependent fixed point of 
(|J) and it is extrapolated to the thermodynamic limit. 
Fitting the form D c l 1 l+2 ~ 1/(^ + 1) to the data, the 
extrapolated value is determined as D c = 0.93 ± 0.01, 
based on the standard least square method. Thus for 
< D < 0.93 the system is gapless in all the region of 
< m < 3/2, while for D > 0.93 the energy gap is in- 
duced just at m — 1/2 and the magnetization curve has 
a plateau. 
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FIG. 1. Scaled gap LA versus the single-ion anisotropy 
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FIG. 2. //-dependent fixed point D c l,l+2 is plotted ver- 
sus 1/L to determine D c in the thermodynamic limit. The 
estimated value is D c = 0.93 ± 0.01. 

The phenomenological renormalization group can also 
estimate the exponent v defined as A ~ (D — D c ) v , using 
the L-dependent form 



r L 



VL,L+2 



= log [—j— ]/log[ 



2 n/l ] {L + 2)A' L+2 {D cL , L+2 ) . 

/ log — Hvm ^ — - 



(5) 



where A' L (D) is the derivative of Al(D) with respect to 
D. The result showed a diverging behavior of vl,l+2 with 
increasing L. It implies that A does not have any alge- 
braic form near D c . Thus the phase transition is expected 
to be the Kosterlitz-Thouless(KT)-type which is also 
consistent with the existence of a finite gapless region un- 
der D c . In addition a naive argument restricting us to 
three states S z =3/2, 1/2 and —1/2 (neglecting the state 
S z = —3/2 because of a large magnetic field) at each site, 
leads to a mapping the Hamiltonian ([!]) to a generalized 
anisotropic 5 = 1 model without magnetic field, which 
has the KT phase boundary between the large-D (sin- 
glet) and XY (planar) phases. J(||l4| 

To determine the universality of the phase bound- 
ary D c at m = 1/2, we estimate the central charge 
c in the CFT and the critical exponent r/ defined as 
(S+S-) - (-l)V-" for D < D c . The CFT || pre- 
dicts the asymptotic form of the ground state energy per 
site as 



\e{l,m) 



7T 1 

e(m) - qCVs-^ 



(L -> oo), 



(6) 



where v s is the sound velocity which is the gradient of the 
dispersion curve at the origin. Thus the central charge c 
can be numerically determined by estimating the gradi- 
ent of the plots of E(L,M)/L versus 1/L 2 and v s . v s is 
estimated by the form jll| 

v s = ^(E kl (L, M) - E(L, M)) + O(-l), (7) 

where k\ = 2n/L is the smallest nonzero wave vector 
for L and (L, M) is the lowest level in the subspace 
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specified by M and k\. The calculated c for D < D c at 
rn = 1/2 is shown in Fig. [| At the boundary D c (= 0.93) 
our estimation gives c = 1.03 ±0.06 and other points also 
have comparable errors. Thus we reasonably conclude 
c = 1 for D < D c . 



Finally using the method in Ref. gfuj, we present the 
ground-state magnetization curve in the thermodynamic 
limit for several values of D; D —0, 1, 2 and 3. For 
D = the system is isotropic and gapless for > m < 
3/2. For other cases, it has the gap at m = 1/2 and the 
magnetization plateau appears. 
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FIG. 3. Estimated central charge c and exponent 77 
for D < D c . At D = D c (=0.93) our estimation gives 
c = 1.03 ± 0.06 and r) = 0.26 ± 0.01. We conclude c = 1 
for D < D c and 77 = 1/4 at D = D c . 

Using another prediction of the CFT A± ~ 
irv s ri/L (L — > 00), the exponent 77 can be estimated 
by the form O] 




L,L+2 



(D) = 



1 + log 



•Ax 



2(D)- 



ATTS) 



/log 



r A^(-P)A' L + 2 (_D) - 
LA £ (r>)Ai +3 (13). 



I3 l .l +2 (D) ~ (D-D cL , L+2 ) 
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The calculated 77 is shown in Fig. ^[ Our estimation 
77 = 0.26 ± 0.01 at D = 0.93 suggests 77 = 1/4 just at the 
phase boundary. In addition the estimated 77 gradually 
decreases with decreasing D. Thus the analysis on 77 also 
supports the KT transition. 

The critical behavior for D > D c can be tested by the 
Roomany-Wyld approximation for the Callen-Symanzik 
/3-function 



FIG. 4. E(L,M+1)-E(L,M) and E(L,M)-E(L,M-1) 
plotted versus 1/L with fixed m for D — 3. The dashed curves 
are guides to the eye. The extrapolated points for m — 1/2— 
and m = 1 /2+ corresponds to the resutls of the Shanks' trans- 
formation H- = 4.17 and H+ = 5.19, respectively. 

Since the system is gapless except for 777, = 1/2, H + 
and 7J_ of correspond to each other and the common 
value gives the magnetic field H for given m in the ther- 
modynamic limit. The size correction of (|^) is predicted 
to decay as ~ 0(1/ L), by the CFT. Thus we can estimate 
H for given to, using the extrapolation form 



E(L, M + 1) - E(L, M) ~ H + 0(1/ L) 
E(L, M) — E(L. M - 1) ~ H + 0(1/L) 



(11) 



(9) 



When the gap behaves like A ~ exp(— a/(D — D c ) a ), the 
function (S) has the form 



(L-kx>), (10) 



in the thermodynamic limit. Fitting the form ( |l0| ) to 
the calculated function (||) for each L, a is estimated as 
follows: cr 8! io = 0.46 ± 0.06, a w ,i 2 = 0.52 ± 0.05 and 
Ci2,i4 = 0.56 ± 0.06. The results are also consistent with 
the standard KT transition (a = 1/2). Therefore we 
conclude the critical behavior near D c for m = 1/2 is 
characterized by the universality class of the KT transi- 
tion. 



with fixed m. For D = 3.0 the left hand sides of the 
form @ calculated for m = 0,1/4,1/2,3/4,1 and 5/4 
are plotted versus 1/L in Fig. |^. It shows that the form 
( |Tl| ) is valid except for m — 1/2 and the two extrapolated 
values of H (the one is extrapolated from E(L, M + 1) — 
E(L, M) and the other is from E(L, M) - E(L, M - 1)) 
correspond to each other well. Thus we take the mean 
value of the two for the magnetic field for each m. Only 
for m = 1/2 H + and H- are obviously different and 
the size correction decays faster than 1/L, as shown in 
Fig. |J, because the system has a gap. Then we estimate 
H + and H- by the Shanks' transformation |T^] P' n — 
(P n _iP n+ i - P^)/(P n _! + P„ + i - 2P„) for a sequence 
{P„}. Applying it twice to E(L, M + 1) - E(L, M) and 
E(L, M) - E(L, M - 1) respectively, for L = 6,8, 10, 12 
and 14, results in H+ = 5.19 ± 0.07 and = 4.17 ± 
0.07, which are indicated as the extrapolated points in 
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Fig. [|. The extrapolated value H for other values of 
D can be estimated in the same way. Only for D = 
H + and H— correspond even at m = 1/2. The ground- 
state magnetization curve in the thermodynamic limit is 
given by all the extrapolated values of H for each m. 
We present the results for D =0, 1, 2 and 3 in Fig. ||, 
where we also used the values of H for m =1/3, 2/3, 5/6, 
7/6 and 4/3 which are estimated by the same method as 
mentioned above. The curve has a plateau at m = 1/2 
(H- < H < H + ) for D =1, 2 and 3, in contrast to 
the case of D = which does not have any nontrivial 
behaviors. 









/ J 
* / 

/ , 


/ J* 
f J* 


• / ,+- 

/ / >■' / 

* / *' y 
4 * # y +- 


— • D=0 


/ / / / 


— * D=1 




-~ + D=2 




— D=3 



0.0 5.0 10.0 

H 

FIG. 5. Ground-state magnetization curves in the ther- 
modynamic limit for D = 0, 1, 2 and 3. 

Among those curves in Fig. |5| D = 1 is the most im- 
portant in terms of experiments to detect the plateau, 
because D ~ J might be realized in some real materials. 
The candidates of the quasi-lD S — 3/2 antiferromag- 
net are CsVCl 3 @ and AgCrP 2 S 6 ||. In particular for 
AgCrP 2 S6 a large anisotropic effect was observed in the 
magnetization measurement in low helds. Higher-field 
measurements of those materials would be interesting. 
Note that for D > D c the ground state is gapless for 
H < H- and H > H + , while massive for < H < + . 
In the quasi-lD systems, some canted Neel orders occur 
in the ID gapless phase, due to interchain interactions. 
Thus a re-entrant transition might be observed in the 
magnetization measurement; with increasing H the Neel 
order disapears at H_ and appears again at H + at suffi- 
ciently low temperatures. 

In summary the finite cluster calculation and size scal- 
ing study showed that the anisotropic 5 = 3/2 has the 
magnetization plateau at m = 1/2 for D > D c = 0.93 
and the phase transition with respect to D belongs to the 
same universality class as the Kosterlitz-Thouless transi- 
tion. 
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